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The high frequency end of the relaxation 
spectrum for polymer molecules involves 
the rotation of the segmental bonds. This 
fast relaxation process, however, cannot 
take place easily in the condensed state 
crowded by the densely packed conformers, 
necessitating the slower cooperatively syn- 
chronous relaxation. As the temperature is 
lowered, the domain of cooperativity grows 
towards the infinite size at the Kauzmann 
zero entropy temperature, though actually 
the system deviates from the equilibrium as 
the glass transition intervenes typically at 
50 K above that temperature. The excess 
enthalpy and entropy drop faster than pre- 
dicted by the rotational isomeric states 
which would reach zero only at K. The 
real ACp is greater than that of the RIS 
value. The actual volume in excess of the 
crystalline lattice volume, however, points 
towards zero at K. Thus, a polymer with 
higher Tg typically exhibits a lower density 



and modulus in the glassy state. Since the 
configuration al entropy associated with the 
free volume is proportional to the logarithm 
of the latter, the Kauzmann temperature 
can be scaled by InM, where M is the al- 
gebraic average of the conformer molecular 
weight. The temperature dependence of the 
most dominant, i.e., the largest equilibrium 
domain size will result in the Adam-Gibbs 
and Vogel equations for the characteristic 
relaxation time. The cooperative domain 
distribution leads to the relaxation spectrum 
that follows a power law. The relationship 
between the characteristic relaxation time 
and the rate of physical aging is derived. 
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1. Introduction 



Is it the relaxation or retardation spectrum that would 
better reflect molecular motions? 

There are two clearly different types of viscoelastic 
experiments. One is to control the strain (deformation) 
and the other is to control the stress. Data of the first 
type include the relaxation modulus G(t), the secant 
modulus under the constant rate of strain, the dynamic 
modulus G *(w), and the dynamic viscosity 17 -^(co ). The 
relaxation spectrum describes the intensity of each 
mode of relaxation with a unique rate of dissipating the 
energy. It is considered to be proportional to the number 
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of relaxing units. The strain energy is closely related to 
the change of molecular configuration in polymers. 

The scheme of relating the relaxation spectrum to 
molecular relaxation seems justified by the apparent 
success of the Takayanagi model [1] for predicting the 
dynamic moduli for heterogeneous structures, such as 
composites, block copolymers and semicrystalline 
polymers, from the moduli of individual components 
weighted according to their volume fractions. 

The second type, or the stress-controlled type, in- 
cludes the creep compliance J(t), the secant compliance 
under the constant rate of stress increase, the dynamic 
compliance /*(w), and the dynamic fluidity (f^^ico). 
Unlike the strain energy, a new approach will have to be 
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devised which will relate the stress energy to molecular 
strains, and this has not been done. More assumptions 
are needed to interpret a retardation spectrum in molec- 
ular terms than a relaxation spectrum. However, these 
data can be converted to the first type mathematically. 

Because the dynamic compliance /*(w) is simply the 
reciprocal of the dynamic modulus G*(w) in one and 
the same experiment, /*(w) can be included in the first 
group, provided the data are analyzed in terms of the 
relaxation, not retardation, spectrum, and its signifi- 
cance can be analyzed in molecular terms. 

In linear viscoelasticity, the creep compliance J{t) is 
not simply the reciprocal of the relaxation modulus 
G{t), but it can be exactly related to G{t) by the equa- 
tion below from the fact that the complex dynamic mod- 
ulus G* = G' + /G" is the reciprocal of the complex 
dynamic compliance /* = T-iJ" : 



t 

I 



G{t-x)J{x)dx = t 



(1) 



together with three other similar convolution integrals. 
If the relaxation modulus can in fact be represented by 
the sum of G/ exp(-?/T/) terms, then it can be seen that 
the creep compliance cannot be in the form of the sum 
of the Ji(l-cxp(-t/Ti)) terms. All models utilizing any 
numbers of springs and dashpots fail to satisfy Eq. (1). 
Springs and dashpots are not suitable in describing lin- 
ear viscoelastic behavior of a given material. If the rate 
process is assumed for the decay of strain energy under 
the constant macroscopic strain, the relaxation modulus 
can be the sum of G/exp(-?/T/) terms, and all other 
viscoelastic functions must be consistent with this as- 
sumption. In such a case the relaxation spectrum, but 
not the retardation spectrum, would reflect the struc- 
ture-property relationship. 

The relaxation spectrum //(In) is approximately: 



//(Inr) = - 



dG(0 
dlnf 



T= t . 



(2) 



It will be shown later that the intensity of relaxation 
//(Inr) is proportional to the number of participating 
conformers, and //(Inr) = G"(lnw), particularly when 
ldlogG7dlogwl«L 

The relaxation spectrum may be "normalized" by the 
absolute value of the modulus, and 



HQnr) ^ G'\(o) ^ r{(o) 



This quantity is approximately equal to the slope of the 
double logarithmic relaxation modulus curve when 
G*(w) = G(t), such as in the solid-like state. The dielec- 
tric loss factor s"(a)) is similar to J"(co), and it too can 
be a basis of molecular interpretation. 



2. The Conformer 

Theodorou and Suter [2] have demonstrated with 
computer simulation that an overwhelming proportion 
of elastic deformation in amorphous polymers is taken 
up by partial rotation of the main chain bonds. This is 
not a surprising conclusion, and it can be carried further 
by suggesting that the irreversible bond rotation from 
one stable configuration to another in the presence of a 
force field is the basic mode of relaxation. The rotation 
of one single bond from a gauche to a trans configura- 
tion is the fastest mode of relaxation involving the 
smallest unit, called conformer. In a vinyl polymer, a 
monomeric repeat unit consists of at least two conform- 
ers. 

When a conformer changes the configuration, it must 
pass an intermediate stage that has a much higher free 
energy level, amounting to 12.6 kJ/mol to 15. 1 kJ/mol of 
conformer. If molecular weight of the conformer is 50 
and the density is 1 g/cm^, this means 3 X 10^ ergs, 
which is 10 times greater than the strain energy spent at 
the yield for a typical glassy polymer. Since the me- 
chanically imposed strain energy is so small as com- 
pared to this barrier energy, the probability for a suc- 
cessful rotation, i.e., relaxation, can be considered as 
controlled by the probability of these thermally excited 
conformers reaching this energy barrier. It becomes the 
bottle neck for the irreversible transition in configura- 
tion, and its relative population determines the rate 
constant for the process. The relative population is 
proportional to the exponential function of negative en- 
ergy difference over kT, according to the Maxwell- 
Boltzmann equipartition principle, so that the rate of 
transition is obtained from the equation: 



/7 0cexp 



kTj 



(4) 



(3) 



where p is the probability that the gauche configuration 
will spontaneously transit to the trans state, and Afx is 
the highest free energy level the bond angle must pass, 
relative to the stable state in equilibrium, k is 
Boltzmann's constant and Tis the absolute temperature. 
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As is well known, this is the classical theory of rate 
processes [3], and the rate constant is proportional to the 
probability p . The time constant for this transition is the 
relaxation time for a conformer and it is proportional to 
/7"\ or exp(A^t/^r). 

This relaxation time is of the order of 10"^^ s at room 
temperature [4], and with A^tof 12.6-15.1 kJ, it is 
estimated to be 10"^^ s at 500 °C. This, then, comprises 
the shortest time end of the vast spectrum of many 
modes of relaxation a polymer molecule is capable of 
undergoing. The dominant "characteristic relaxation 
time" of polymers is many orders of magnitude slower 
even at 50 K above Tg. In the dense liquid with crowded 
conformers, relaxation requires a simultaneous rotation 
among several neighboring molecules, and this can in- 
crease the relaxation time by many orders of magnitude. 
Such cooperativity is not unique to polymers, but gen- 
eral to most glass-forming liquids. 

Polymers are unique, on the other hand, in their capa- 
bility for exhibiting an extremely slow relaxation pro- 
cess well above their glass transition temperature, and 
this behavior too has been described by the same word 
"cooperativity." The context for this latter usage is dif- 
ferent from the former, as the latter involves a large 
overall configurational change of strings of conformers 
which make up a polymer chain. Their individual mo- 
tions add up to a net movement of a large segment of a 
chain. Coordinated motions among strings or "strands" 
of a polymer molecule can be slower by orders of mag- 
nitude than the cooperative relaxation among the indi- 
vidual conformers described above, though all relax- 
ation modes in the latter context are coupled to ones in 
the former context. The former governs the time-tem- 
perature shift of the entire spectrum including the latter. 
The elastic modulus for the latter arises from a decrease 
in molecular conformational entropy accompanying the 
deformation, a typical rubber, whereas the elastic mod- 
ulus for the former arises from the increase in internal 
energy, a typical solid. 



of the relaxation process as we understand it, begins to 
be lost above 7*. We arrived at the value of 500 °C for 
T^ after analyzing data for more than a hundred poly- 
mers [1]. 

If conformers were able to rotate without being inter- 
fered with by their neighbors, the activation energy bar- 
rier would remain constant at A^t, ca. 12.6 kJ, and the 
temperature dependence of the relaxation time r (which 
is very short) would follow the Arrhenius relationship: 



Inr = Inr* + 



kT 






(5) 



where r* denotes the relaxation time at 7*. 

Because of the neighboring interference, however, the 
relaxation time grows at a faster rate with temperature 
change than shown by Eq. (5). The degree of coopera- 
tivity is characterized by defining the domains of coop- 
erativity, in which all conformers must relax simulta- 
neously. The domain size, z , then becomes the exponent 
for the probability for simultaneous relaxation of other- 
wise independent conformers. We now obtain, instead of 
Eq. (5): 



Inr = Inr* + 



kT 






(6) 



Equation 6 is no longer an Arrhenius equation because 
z changes with the temperature. The concept of the 
temperature dependent domain size, z, is the center 
piece of the theory by Adam and Gibbs [5], although 
their further treatment of it seems to imply the in- 
tramolecular, rather than intermolecular, cooperativity. 
The property of z is such that z= 1 at 7*, and z^^ 
at To, where the domain size has grown so large that all 
conformers are locked together. Such a temperature de- 
pendence of z can be met with the equation: 



z = 



T-To 



T^-To 



(7) 



3. Intermolecular Cooperativity 

Our discussion will be limited to the "cooperativity" 
in the first context, i.e., the relaxation process in which 
a number of conformers relax simultaneously. At a high 
enough temperature, the volume is expanded to allow 
most conformers to relax individually. In the actual 
polymers this may be well above the degradation tem- 
perature, but the extrapolated value will be an important 
parameter. This temperature is defined as T^. At this 
temperature, the frequency and the magnitude of oscil- 
lation are such that a change to a new bond angle can be 
reversed in negligible time, and the theoretical meaning 



A real liquid vitrifies at about 50 K above To, but 
amorphous glassy polymers are far from the state of 
zero configurational entropy [6]. Extrapolation from 
above T^ to 50 K below T^ will hardly increase the 
molecular order, and it is safe to say that configurational 
entropy is unlikely to extrapolate to at Tq. The config- 
urational entropy ^c should be estimated from the statis- 
tics of the rotational isomeric state, and this entropy 
would reach at K. So, what kind of entropy is it that 
Kauzmann [7] observed heading toward at TqI To 
clarify this question, one needs to examine how the 
configurational entropy per mole of conformer changes 
as the domain size grows. 
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The configurational entropy ^c is calculated from the 
concentration dependence of the gauche and trans con- 
figurations, where the concentration cj) for the gauche 
state depends on the difference in the energy level from 
the trans state, AF: 



cj) = 2Qxp(-Ar/kT)/[l + 2 exp(-Ar/^r)] 



(8) 



The factor 2 is placed in the formula because there 
are three possible angles that a conformer can assume 
per bond. In some polymers, the third position is forbid- 
den by the geometrical interference; in such a case this 
factor 2 is replaced by unity. The entropy of mixing 
between the trans and gauche configurations is obtained 
through an approximation: 



Sc = -6k[cl)\ncl) + (1 - (/))ln(l - cj))] 



(9) 



The factor 6 is placed in the equation because each of 
the two conformers (trans or gauche) has a choice of 
three bond angles. 

Now, in a domain that consists of z conformers, all 
configurational changes must occur together, and that 
means the degree of freedom inside the domain is the 
same as that of one free conformer. The configurational 
entropy ^c is reduced to S^ = SJz because of the cooper- 
ative restrictions. In other words, the entropy 5c is 
reduced at all temperatures but still headed toward zero 
at K, while Sz = SJz decreases rapidly toward zero at 



Tq. This is illustrated in Fig. I, showing ^c, SJS, SJ6, 
and Sz- 

Although the celebrated theory of Gibbs and Di- 
Marzio [8] is a theory that leads to zero entropy at To, 
it is different from the cooperative domain model pre- 
sented here. Gibbs and DiMarzio modified Flory's lat- 
tice model, but by including holes that are temperature 
dependent, the conformational entropy is made to reach 
zero at Tq. Still, their theory requires that polymer 
molecules be sufficiently ordered near To, whereas with 
the model presented here the configurational entropy is 
reduced by virtue of the increasing cell size of the 
microstates as the domain grows at lower temperatures. 
Since the glass transition is not unique to polymers, a 
model must be easily adaptable to nonpolymers. The 
Gibbs free energy can be calculated from the tempera- 
ture dependence of Sz , and from it the enthalpy Hz can 
be, also. Figure 2 is an illustration of the temperature 
dependence of TSz, Hz, and the (negative) Gibbs free 
energy; Hz < TSz for Tq<T < T*. The plot of TSz is 
nearly a straight line starting at zero at To, which is 
useful in scaling among polymers with different values 
of To. 

The specific heat, ACp, i.e., the difference between 
the liquid and the solid (glass) state, is of interest, and 
is shown in Fig. 3, together with that which corresponds 
to the 5c. Both values are per mole of conformers, and 
they must be divided by the molecular weight of con- 
former to obtain the value per g. The peak for dHJdT 
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Fig. 1. The confomiational entropy per conformers, per domains of size z = 3,z = 6, and z tliat depends 
on temperature (the dotted line). Tq = 370 K is chosen. (Note that the units of entropy are in moles of 
conformer.) 
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Fig. 2. Enthalpy Hz, the product TSz, and the Gibbs free energy (negative value), all per mole of 
conformers, calculated from the 5z = SJz{T) shown in Fig. 1. 
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Fig. 3. The specific heat ACp per gram has been calculated since the molecular weight of the 
conformer can be calculated from Tq. This result agrees with the data for polycarbonate. 



appears as the consequence of the increased order at 
lower temperature, and the greater the energy difference 
AT" between the trans and the gauche configuration, the 
greater this peak appearing at a higher temperature, but 
this should not be confused with the glass transition. 

ACp depends on temperature more strongly than ACp 
cc 1/r, a prediction by Scherer [9]. It will be shown later 



that the specific heat per mole is proportional to the 
ratio T^/(T^-To), and to InM, where M is the molecular 
weight of a conformer. ACp per mole of conformer is 
thus greater for a polymer with a larger conformer, 
which has a higher Tg, but Cp per gram is smaller for the 
same with the larger conformer, following the formula 
(lnM)/M, with higher Tg. 
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The relationship between S^ and Sc can also be found 
parallel to nonpolymeric liquids that are capable of ex- 
hibiting internal configurational changes, known as 
Kneser liquids [10]. Domains of cooperativity among 
the molecules can be considered, within which the 
molecular relaxation can only occur together simulta- 
neously. By the same token, S^ = SJz decreases toward 
zero at Tq. Herzfeld and Litovitz [11] find the absorption 
of sound waves in these Kneser liquids to be many times 
greater than the estimated "classical" values. A sound 
wave consists of the adiabatic compression and expan- 
sion, and dissipation of free energy occurs when heat 
conduction irreversibly equalizes the temperature differ- 
ences between the crest and the trough. 

The "fragile liquids," according to the well-known 
classification of liquids by Angell [12], are also those 
which can exhibit internal degrees of freedom when 
free from the neighboring interference. In almost all 
"fragile" cases, \dz/dT\ is large, i.e.. To is closer to T*, 
and r*/(r*-7'o) is large, and the relaxation time follows 
the Vogel equation. The other class of liquids termed 
"strong" by Angell is presumably that for which the 
domain size either does not grow or grows little over the 
temperature span, and typically exhibit Arrhenius 
behavior. A "strong" glass has less free volume than a 
"fragile" one. 

Substituting the temperature dependent z, (Eq. 7), 
into the Adam-Gibbs equation [Eq. (6)] will immedi- 
ately give the Vogel equation [13]: 



lnT= Inr* + 



A/i^ 



A/i^ 



k(T-To) k(T^-To) 



(10) 



paraphenylene linkage has about the same value. The 
ether linkage, on the other hand, has only 4.8 kJ, so an 
adjustment should be made when estimating T^ from 
polymers that include it in the repeat unit, as will be 
discussed later. 

The Vogel equation, Eq. (10), is exactly the same as 
the Williams-Landel-Ferry equation [14] with appropri- 
ate substitution of the parameters, and also the same as 
the Doolittle free volume Eq. (15) when the thermal 
expansion coefficient of the fractional free volume (an 
empirical parameter different from the van der Waals 
free volume) is equated to A:/A/i* in the Vogel equation. 
This will be shown later. For the "universal" WLF 
equation, C2^ = Tg - To ^ 50 K, and d^ = 17.4 [16]. The 
latter is related to the universal A/x*, while the former 
is an indication that a typical glass transition occurs at 
50 K above Tq. All these equations pertain to the condi- 
tion for thermodynamic equilibrium, i.e., above T^. 

To summarize the relationship among these quanti- 
ties and limiting temperatures, the following equation 
may be written: 



z = 



TS, 



r* T-To 



(12) 



Substitution of Eq. (12) into the Adam-Gibbs equa- 
tion, Eq. (6), results in another, more widely known 
form of the Adam-Gibbs formula: 






-^^ (13) 



where the relationship between and A/x and A/x* are 
related by the formula: 



A/i = A/i^ 



7^* 



T^-To ' 



(11) 



The ratio r*/(r*-ro) above has appeared already, 
and is an important empirical material parameter. 
Recalling that the product TS^ in Fig. 2 is a straight line 
which reaches zero at To, the ratio {TSc/TS,} = T*/ 
(X^-T^o) can be scaled for all materials with different To 
values. The physical significance of the ratio T*/ 
(X^-T^o) is the conversion of the configurational proba- 
bility from without intermolecular cooperativity to with 
intermolecular cooperativity. All thermodynamic ex- 
tensive quantities are scaled by this ratio, including ACp 
and free energy such as Afi . The true rotational barrier 
of a bond is A/x* when it is free from the neighbor 
interference, hence A/x* is a universal constant for con- 
formers that involve a C-C bond rotation with the value 
of a little more than 1 2.6 kJ/mol. Conveniently, Afi * for 



4. Free Volume and Entropy 

Equation (10) is a special Vogel equation for the shift 
factor from T*. If T^ is chosen instead of T* as the 
reference temperature, Eq. (10) is changed to: 



lnT= Inxe + 



Ajji"^ 



A/i'^ 



'^k(T-To) k(T,-To) 



(14) 



The above equation can be written in the form of the 
Doolittle' s free volume equation. 



In T = In Tg + 



1 



1 



af(T-To) af(Tg-To) 



(15) 



where a/is the thermal expansion coefficient of the free 
volume fraction / = (^-^0)/^, where Vo is the 
"occupied volume" V at Tq. 
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We have already pointed out that a^is really a parame- 
ter, A:/A/i*, arising from the energy barrier for bond 
rotation, A/x*. The smaller the A/x*, e.g., for the ether 
linkage, the greater the empirical a/, irrespective of the 
thermal coefficient for the van der Waals (real) free 
volume, Aa, the difference in a between the liquid and 
the solid (a glass or the crystal). Whereas the Doolittle 
free volume vanishes at Tq like 5^, the "real" free volume 
vanishes at K like S^. In order to be consistent, the 
following relationship should hold: 



Aa 



T^-To 

To 



(16) 



and here again we find the familiar ratio. The empirical 
value for a/ ~ R/Afx'^ obtained from viscoelastic data is 
about 6 X 10^ K"^ if A/x* is a Httle more than 12.6 kJ, 
whereas the values of Aa and a/ for many amorphous 
polymers vary around 4 and 6X10^ K"\ respectively 
[16]. The ratio Aa/a^ ought to be 0.5 if T^ = 150 °C, and 
0.85 for Tg = - 1 10 °C. For polyethylene, taking a for the 
crystalline lattice for a of the solid, the value for Aa = 
ai-as is 5.7 X 10^ K"^ [17]. Furthermore, if one as- 
sumes A/i* to be 16.7 kJ, Doolittle's /g of 0.025 is 
obtained. 

The "real" van der Waals free volume with the ther- 
mal expansion coefficient of ca. 6X10"^ K"' will attain 
50 % of the van der Waals volume at T*. Thus a true 
excess volume is comparable to the true occupied vol- 
ume of a conformer in the state free from the neighbors. 
A greater amount of free volume will be trapped in the 
glassy state if the Tg is higher, as has been found for 
crosslinked polymers with different degrees of 
crossHnking density [18]. Such a trapped free volume is 
a genuine unoccupied site that can be observed by the 
positron annihilation analysis. As a rule, the glassy mod- 
ulus of a higher Tg polymer tends to exhibit a lower 
value, which is further evidence for the less dense glassy 
state for the high Tg polymers [19]. 

Now, the configurational entropy Sf of the molar free 
volume Vf can be considered to be equal to the configu- 
rational entropy of the ideal gas with the molar volume 
of Vf, i.e.. 



Sf=-k\nVf 



^f 



/■ 



(17) 



The absolute value of the molar free volume scales 
with molecular weight of conformer M divided by the 
density. Since the value of the excess entropy at T* is 
the same for all polymers, the product (r* - ro)ln Vf is 
constant for all polymers. When polymers with different 
conformer sizes are compared, their values of T*- To 
are inversely proportional to ln(M/p), where M is the 



molecular weight of the conformer and is the density 
extrapolated to T*, or 



J-* 



InM 



T^-To 



(18) 



Here again, we find the same ratio T^I{T^ - To) for 
scaling the extensive quantity. 

Ignoring the density variations among the amorphous 
polymers, (in comparison to the molecular weight varia- 
tions among different conformers), the constant 
C = (r* - To) In M has been evaluated through the intro- 
duction of a virtual polymer with molecular weight of 
Mo, such that its To is at K, i.e., z = 1 and S^ = Sc for 
all temperatures from K to T*. We obtain the formula: 



(r* - To) InM = r*lnMo = C . 



(19) 



This relationship is shown in Fig. 4. 

The empirically derived value for C is 1750 for many 
polymers, which means Mo = 9.6, too small for a 
methylene unit and not even as large as one carbon 
atom. For the methylene unit as a conformer, T^ of 
- 1 10 °C is obtained from Eq. (19). Because this is near 
the y transition temperature of branched (low density) 
polyethylene, it is often regarded as the same jS-type 
relaxation. Their temperature coefficients are substan- 
tially different [20], suggesting that the glass transition 
for the linear polyethylene and the y transition for the 
branched polyethylene overlap in the kHz range. The y 
transition for branched polyethylene is definitely a local 
relaxation process involving a fixed number of 
intramolecular cooperative conformers, exhibiting a 
constant activation energy, just as in the case of the jS 
transition in amorphous homopolymers, as will be 
discussed subsequently. 

The above discussion, including tables of Tg for a 
hundred polymers, has been included in Chap. 2 of Ref. 
[20], but it has been further elaborated here to make the 
present discussion complete. Equation (19) is useful in 
predicting the To, and therefore Tg --^ To + 50 K from the 
chemical structure of an amorphous polymer. The gen- 
eral rules for estimating the conformer size from chem- 
ical structure is elaborated below: 

1 . A conformer with a branch. 

(a) When another conformer is attached as a 
branch to a conformer in the main chain, the 
new combination attains the combined size as 
one conformer. When a methyl group is at- 
tached to a carbon in the main chain, this 
combined conformer consists of CHCH3 and 
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Fig. 4. The glass transition temperature, 7^ = Tq + 50 °C, vs segmental or 
conformer molecular weight, according to Eq. (19). 



its molecular weight is 28. However, 
polypropylene has a methylene unit as the other 
conformer along the main chain. Averaging the 
two, (28 + 14)/2 = 21 is obtained, and To of 
-75 °C, or Tg of -25°C is obtained from Eq. 
(19). 
(b) When the branch contains more than one con- 
former so the branch itself can undergo its own 
configurational relaxation, the number of con- 
formers per monomeric unit increases and the 
average conformer size decreases again. For 
polyvinyl acetate, the monomeric unit consists 
of three conformers instead of two. The average 
conformer molecular weight M is 86/3 = 28.7, 
and from Eq. (19), Tg of 29 °C is obtained. 
Whether a further addition of conformers will 
raise or reduce the Tg depends on the size of a 
conformer to be added. Simply, if it is larger 
than the present average size, the Tg will rise, 
but if smaller T, will fall. 



2. The end group in a polymer chain has a higher 
enthalpy and entropy than a similar conformer inside 
the chain. The end group lowers Tg as would the 
solvent molecules. The contribution of the chain 
ends to the overall enthalpy can be estimated by 
mixing them with those conformers along the chain, 
and the average Tg can be obtained. 

3. When compatible polymers are mixed, a value of Tg 
between the two individual values is found. One 
approach for estimating the T^ from the individual 
concentrations of the two polymers is to invoke the 
principle of isentropic state at Tg = To + 50 K. The 
entropy ^g for cooperative relaxation determines the 
domain size z at Tg, so the same value for ^g means 
the same z value at T^ for all polymers and their 
mixtures. The entropy ^g can be derived by the 
formula: 






(T,-To) 



(20) 
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but Tg - To is 50 K, and dS/dT at T^ is A Cp per mole of 
conformer, so we have: 



„ 5QACp 



(21) 



and A Cp per mole is proportional to T'^/(T'^-To) or to 
InM, hence 



S ^ 



InM 



(22) 



is obtained. A mixture containing x mole fraction of 
conformer A and 1 - x mole fraction of conformer B 
will exhibit the average Tg, 



InM InMA ^, , InMg 
-7^ = x— — + (1-x)— ^ 

J g J gA J gB 



(23) 



which reduces to the Fox-Rory formula [21] when the 
variations in InM are negligible. 



5. The Relaxation Spectrum 

The discussion so far has been with a single value of 
z that depended on the temperature. There no doubt is, 
however, a dynamic distribution of domain sizes at any 
temperature. The domain size as discussed up to this 
point is the most dominant size, which from this point 
on will be called the "characteristic size" Zc, and the 
characteristic relaxation time r^ corresponds to this Zc Tc 
is the relaxation time in the Kohlrausch- Williams -Watts 
formula [22], exp{-(r/T)^}. In the high temperature 
range, e.g., more than 100 K above Tg, Zc is small and 
the spectrum is narrow. At moderate temperatures, e.g., 
50 ~ 100 K above Tg, the spectrum broadens and the 
Cole-Cole parameter ^ [23] is typically 0.6 to 0.9. At 
even lower temperatures (but in equilibrium above Tg), 
the spectrum broadens particularly on the high fre- 
quency side, signifying the increase in domains that are 
smaller than Zc, primarily because Zc has become larger. 
Cole-Davidson [24], KWW, and Havriliak-Negami 
[25] formulas are all for "nonsymmetrical" distributions 
around the characteristic relaxation time for the loss 
peak. (Empirically, Ngai's [26] stretched exponential 
function is equivalent to the KWW formula.) When 
plotted in a doubly logarithmic scale, all three becomes 
a straight line in the high frequency side, while only 
H-N gives a straight line on the low frequency side with 
an additional adjustable parameter. Polymers exhibit an 
additional weak but broad extension of the relaxation 



spectrum on the low frequency side that follows a 
power law, with a value between 0.5 and 0.8 for the 
exponent, which is due to the slow rearrangement of 
molecular chains known as the Rouse-Zimm type distri- 
bution, as briefly mentioned in Introduction; the spec- 
trum being discussed here arises from density fluctua- 
tions causing a dynamic distribution in domain size. 

Ngai et al. [26] assume that the relaxation spectrum 
is homogeneously broadened by coupling to low energy 
excitations. Chamberlin et al. [27], on the other hand, 
consider a heterogeneous distribution of exponentially 
relaxing domains, similar to the domain size distribu- 
tion discussed here. Chamberlin's model assumes a 
Gaussian distribution of domains, which results in re- 
ducing the slope of log e" vs log ft at very high fre- 
quencies. The domains we consider are from the hetero- 
geneity bom of density fluctuations. Those domains 
with sizes that are far different from the equilibrium 
size, Zc, have higher free energy, and are more scarce at 
both extremes. We apply the equipartition principle to 
estimate the distribution of the number n^ of conformers 
that belong to domains of size z : 






exp 



A^ 
" kT 



(Zc-z) 



(24) 



where ric is the number of conformers in a domain of 
size Zc, and A ^is the Gibbs free energy increase when 
the domain size is disturbed from Zc by A z = 1 . 

From the Adam-Gibbs formula, Eq. (6), the formula 
for the perturbation in z from the most dominant size Zc 
can be obtained: 



— = exp 

Tc 



[-^(zc-z)] 



(25) 



and 



Tz IS 



where Tc is the characteristic relaxation time, 
the relaxation time corresponding to size z . 

By combining Eqs. (24) and (25), we obtain the for- 
mula for the relative number of conformers with relax- 
ation time T^, which is the relaxation spectrum H(\nr): 



( \ — 
//(.nx) = (^)- 



c/ 



(26) 



which is a "power law" type formula, i.e., the plot of log 
G"or log e" vs log il is a straight line with the slope 
-A '^/A/i. A '^is the amount of free energy increase 
when the domain size is decreased by 1 , and its value 
should be close to the difference between the gauche 
and the trans configurations, typically 4.8 kJ, so the 
slope is around 1/3. Dielectric data on polycarbonate 
obtained in the temperature range from 155 °C to 
200 °C is shown in Fig. 5. This curve was constructed 
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-2-10 1 2 3 4 5 6 

log freq Hz 

Fig. 5. The master curve of dielectric loss spectrum for polycarbonate, obtained by shifting isothermal 
data. 



by shifting the individual isotherms so the loss peak is 
located at the origin. The slope on the high frequency 
side is 0.3, which would be obtained if A t/f is 4.8 kJ and 
A^t 14.6 kJ. The slope of similar value is observed in 
polyvinyl acetate and a -glucose. The low frequency 
side of the plot in Fig. 5 is steeper, about 0.5, as it is the 
beginning of the Rouse-Zimm spectrum. We believe, 
however, that the Rouse-Zimm model should apply only 
to polymer solutions, and that the relaxation in bulk is 
remote from such a model. A more appropriate model 
for the polymer bulk can be constructed from the distri- 
bution of strands that are dynamically active [20]. Since 
the shorter strands are more numerous, their contribu- 
tion to the elastic energy is greater, whereas the longer 
strands will cause the viscosity to be greater, in propor- 
tion to the strand' length. This results in the relationship 
that the relaxation time increases as the square of the 
strand length, while the modulus intensity is inversely 
proportional to the strand length, and the value of 0.5 is 
obtained for the slope of the spectrum on the log-log 
plot. Another model has been proposed by Tobolsky 
[28] who assumed the normal distribution for the transi- 
tion zone, similar in concept to Chamberlin's model, 
and obtained the slope of 0.5 for styrene-butadiene rub- 
bers of various compositions. Tobolsky notes that stiffer 
chains would affect the elastic constant more severely, 
so that the slope could be greater. A value as high as 
0.75 was found for some systems; this scheme therefore 
can include the Zimm value as well. 

In Fig. 6, the Vogel frequency [= 1/2 ttTc from Eq. 
(10)] is plotted for the polymer with A^t* of 14.6 kJ and 



Tq of 97 °C (dotted line) together with the Arrhenius 
plot by fixing domain sizes in the Adam-Gibbs equation 
for z, Eq. (6). Often, the slope of the Vogel line (multi- 
plied by 2.3/?) is quoted erroneously as the "activation 
energy" with values amounting to 480 kJ or more, ex- 
ceeding the energy for chain scission. The derivative of 
the Vogel equation with respect to MRT is really the 
sum of the correct activation energy A^t * z and the term 
A^t- [dz/d(l//?r)]. The latter is the rate of change of the 
domain size with temperature; this can be a substantial 
quantity. A physically meaningful activation energy is 
the first term z, and this is about 96.2 kJ at 165 °C. At 
147 °C, it is about 12.6 kJ, which is about the right value 
for viscoelastic behavior of glassy polycarbonate. The 
spectrum shown in Fig. 5 is a snap shot at 165 °C, and 
denotes the relative intensity of these Arrhenius lines at 
1000/r~2.28. The intensity for smaller domains de- 
creases with decreasing temperature, at the expense of 
newly created larger domains. Below the glass transi- 
tion, the size Zc becomes fixed, and the characteristic 
relaxation time now follows the Arrhenius line for that 
z, such as a straight hne marked as 3.36, Zc at Tg. The 
shape of the spectrum becomes markedly broader below 
Tg as the activation energy is different for each domain 
size. 

6. The j8 Relaxation 

The conformers within the repeat unit of a chain are 
typically of different sizes. When larger conformers 
become pinned by the neighbors (belonging to other 
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Fig. 6. log frequency from the Vogel equation, Eq. (10), in dotted line. Five straight lines are from 
the Adam-Gibbs formula Eq. (6) with fixed z- At Tg, it is estimated that z = 3.36. 



chains) at the glass transition, there still remains a pos- 
sibility for local relaxation by the smaller conformers in 
the same chain. As a rule the p relaxation requires the 
intramolecular cooperation among adjacent conformers 
in the same chain. The number of cooperating conform- 
ers corresponds to the number of conformers between 
the pinned ones. It is usually the number of conformers 
per repeat unit. The domain size zp for the p relaxation 
does not vary with temperature, and it is equal to the 
ratio of the molecular weight of the repeat unit over the 
average molecular weight of all conformers. This is 
shown in Fig. 7 for polycarbonate. The number of con- 
formers per repeat unit is 3.3, for which T^ of 147 °C is 
obtained from Eq. (19). (The fraction is a result of 
polycarbonate containing ether linkages which requires 
considerably smaller activation energy.) The predicted 
value of Zi3 is thus also 3.3. The activation energy, 
AfjL * * z^ , is 46 kJ from our dielectric data, which sup- 
ports the value of 3.3 for z^ [20]. The formula for the p 
relaxation is given by: 



lnT, = lnT*+^(i 






(27) 



Here A^t* is used as it is simply the energy barrier 
without need to adjust for the density that affects the 
intermolecular cooperativity. The intensity of the p 
relaxation is affected by the temperature and the 



frequency, as it depends on the population of conform- 
ers pinned in the domains. In the high temperature/fre- 
quency region where the a and p lines meet, the p will 
take over because the highest possible frequency is lim- 
ited by the intramolecular cooperation even after the 
intermolecular relaxation has become easy. This is 
shown in Fig. 8 by enlarging the range of temperature/ 
frequency where the a and p peaks merge, for polyvinyl 
acetate. The repeat unit in the polyvinyl acetate 
molecule consists of three conformers, one of which is 
in the branch. In this case, z^ is clearly that of the 
number of conformers in the repeat unit, as A^t*z^ is 
-40.6 kJ and A^t* is 13.8 kJ. 



7. Physical Aging 

It is well known that the glass transition is not a 
thermodynamic transition between two equilibrium 
states, and that a glassy state is a nonequilibrium state 
born of too slow a rate in rearranging the molecules to 
reach the equilibrium state. Physical aging is a process 
for such a glassy state to further shift toward a state of 
lower enthalpy, entropy and free volume. The most 
curious but practically important question is the quanti- 
tative estimation of the characteristic relaxation time Tc 
after the aging period 4 at a temperature T below the 
glass transition temperature. 
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Fig. 7. log frequency for the j8 process which follows an Arrhenius equation, Eq. (27). The dotted line 
is from the Vogel equation for the process as shown in Fig. 6. 




2.75 



Fig. 8. A detail of the log frequency vs l/T for polyvinyl acetate in the frequency-temperature region 
where the a and j8 processes meet. Above the merging point, the process, an intramolecular cooperative 
relaxation process, takes over. 
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Assuming that the aging is a slow process as com- 
pared to the temperature equilibration by thermal con- 
duction, the intensive quantities such as pressure and 
temperature are considered to equilibrate quickly, say in 
seconds, while the domain size z takes a much longer 
time to equilibrate. Consequently, all those extensive 
quantities that depend on z will be slow to equilibrate. 
Thus, in the Adam-Gibbs equation, z and Sz are time 
dependent while the temperature T is not. This makes 
the relaxation time Tc also time dependent. 

If the liquid is vitrified at temperature 7}, the domain 
size thereafter becomes fixed as Zf , which was the 
equilibrium value of z at 7}, i.e.. 



^/ = 



r* - To 



Tf-T,, 



(28) 



Upon further cooling, Zf first remains unchanged but 
will slowly increase with physical aging. 7} is the fictive 
temperature which characterizes such a glassy state. 7} 
is equal to the temperature at the intercept of the volume 
curve for the glassy state and the liquid state. 

Substitution of Zf in Eq. (28) into the first Adam- 
Gibbs equation, Eq. (6), will result in the equation for 
the nonequilibrium state: 



In 7/; r = In 7/; 7^ + 



A^ 



If. 



Ajit^ 



kT Tf-To k(T^-To) 



(29) 



Unlike the Doolittle free volume/, this Xf^r is a physi- 
cally acceptable form of a function in 7} and T. Struik's 
[30] modified free volume is essentially this Z/r- 

In Kovacs' [31] well-known isothermal aging experi- 
ment the sample was initially made to equilibrate at Tini 
above Tg, then suddenly brought to a bath maintained at 
the aging temperature Ta, and the volume was measured 
with aging time 4. Most of his experiments were con- 
ducted near Tg, and the initial relaxation time is in the 
range of (equilibrium) Tg, on the order of tens of sec- 
onds. Not much will happen to Zy^ ^ even after 10 seconds 
has passed, however. There seems to be yet another 
induction period before X/^t begins to change. As we 
shall see, this is because the rate of rearranging domains 
is slower (by more than an order of magnitude) than the 
rate of viscoelastic relaxation. 

Now, in all linear viscoelastic functions, in which t 
and T always appear together as a single reduced variable 
t/r, the relationship shown below is true: 



dlnr 
d\nt 



1 . 



(32) 



The steady state for the aging process is reached when 
the viscoelastic relaxation time T/^t and aging time 4 can 
qualify as r and t in Eq. (32), which can also be written 

as: 



and Tf,r is the characteristic relaxation time Tc with 7} at 
r, and is time-dependent through 7}. This expression has 
essentially been proposed by Hodge in 1986 [29], who 
was able to reproduce A Cp data for polycarbonate sub- 
jected to various thermal histories. This equation obvi- 
ously reduces to the Vogel equation, when 7} is equated 
to r, which is the condition for equilibrium. The Vogel 
equation is good only for the equilibrium condition, 
because the extensive and intensive thermodynamic 
quantities in it cannot be separated. The Doolittle free 
volume equation is also good only for the equilibrium 
condition, as both T and 7} are imbedded inseparably in 
its empirical free volume parameter/. 

Let us now introduce, for convenience, a thermody- 
namic variable X/^t, 



Xfj - 



kT Tf-To 



A)a* 



(30) 



such that Eq. (29) can be written as: 



1,, = ___^1,,. _______ 



(31) 



dr 
dt 



(33) 



At the beginning t^ is zero, while t^t has the initial 
value, so the initial rate of increase for r would have to 
be Qo; 4 cannot be equal to t above, r initially waits 
without changing much until 4 increases enough to 
satisfy Eq. (33). This is the induction period mentioned 
above, and it is to be distinguished from the period of 
temperature equilibration. If a linear viscoelastic exper- 
iments were conducted during the induction period, the 
data will not show the shift according to Eq. (32), i.e., 
one decade of elapsed time resulting in the shift of the 
viscoelastic data by one decade. Rather, the shift would 
be much less than one decade. 

The induction period starts with t^<Tfj, and dlnr/dln/ 
< 1, until the proper ratio of rlt as dictated by X^j is 
arrived, whereupon the ratio T/,r/4 is maintained. For 
example, if the sample is subjected to an initial deep 
quenching, the initial r is very large. Nothing will hap- 
pen at first as the rearrangement of structure is 
extremely slow. After a long wait, t^ has elapsed suffi- 
ciently to the point Eq. (32) is satisfied with 4 as /. At 
this point, X^t, Tfj , 4, and the rate of change of Z/,r all 
become consistent among themselves. The steady state 
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has now been reached, and it will continue indefinitely 
if the aging temperature has been chosen to be below Tq. 
If the aging temperature has instead been chosen to be 
above To, then Xf^r will eventually arrive at an equi- 
librium value where 7}= T^. t/;^ will not increase further 
while 4 continues to increase, and Eq. (32) is again not 
applicable. The steady state condition exists between 
these two extremes, and here the nonequilibrium Adam- 
Gibbs formula can be applied. The kinetics of the steady 
state aging process can be derived by utilizing these 
relationships among X/^, In r. In t, etc.: 



1 


dfa 


1 


dT/,r 


ATfj 
dfa 


-^!. 


dlni/.r 

^ dTfT 


dT/,r 
df. 







(34) 



In deriving Eq. (34), no assumption was made to let 
Tfj~^ be the "reaction rate constant," i.e., we made no 
assumption that the apparent aging rate is the same as 
the rate of configurational relaxation, and in fact it is 
not. The solution of this differential equation is an expo- 
nential integral, and the plot of X/,^ has the appearance 
of a stretched exponential even though a single relax- 
ation time was assumed for the viscoelastic relaxation. 



If a viscoelastic relaxation experiment is started at 
which point it has been aging for 4, G (4) will decrease 
to Git^le in Tfj seconds. Reduction of X/,^ at 4 to \le of 
this value will take much longer, approximately TfjlXfj. 
Xf^T is a small number, ca. 0.1 to 0.01. 

We know the condition at the beginning of steady 
state. The initial value of X/^t requires tedious calcula- 
tions, but it is not far from X^{TIT ), and the ratio of Tfj 
to the aging time 4 is given by the formula: 



^^yI^t 



k(T,-To) 
Afi 



(35) 



Equation (35) is applicable even to the case when the 
initial temperature is below 7g as in the "memory effect" 
experiment and also most of annealing experiments 
where the polymer is heated to below but not too far 
from Tg. The value of the last term is numerically close 
to the Doolittle free volume fraction /g , which is about 
0.025. More detailed approximation yields values rang- 
ing from 0.1 to 0.01, depending on the interplay of the 
initial X and r. Experimental stress relaxation data sup- 
porting r/t^-^fg can be found in Ref. [20]. 

Dielectric data for tan 8 during isothermal aging, as 
measured at 1 kHz, are shown in Fig. 9. The experiment 
followed the procedure that: the sample is kept at 80 °C 
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Fig. 9. Tlie dielectric tan 8 for polystyrene as measured at 1 kHz at 80 °C while isother- 
mally aged. The abscissa is the time starting at each measurement. The "days" in the figure 
means the aging time at 80 °C, and the number in brackets is the log(aging time in 
seconds). 
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while measurements were started from time to time at 
various stages of aging. The time on the abscissa is the 
time spent for each dielectric measurement starting at 
the specified aging time, so the abscissa is comparable 
to -In n . The dielectric tan 8 peak is usually located at 
a higher frequency than for the loss peak s", so the 
difference between the time the peak appears and the 
aging time 4 should be greater than the difference be- 
tween log T and log t^. It is about two decades. The 
dielectric loss peak appears at T^ when measured at 
1 00 Hz, and 1 kHz measurement would make little dif- 
ference, ca. less than 2 °C in the temperature of the loss 
peak. The curve for the longest aging period (120 d) 
suggests that the equilibrium is approached. The relax- 
ation time begins to stay constant while the aging time 
is increased indefinitely. The shortest period of 0.25 d, 
on the other hand, indicates the quenched condition 
which requires some aging time before the relaxation 
time begins to shift. Between these two extremes, the 
ratio in Eq. (32) remains nearly 1 . 

Another caveat for Eq. (32) has not been mentioned, 
and this concerns the distribution of domain sizes, 
which will "stretch" the spectrum. The aging process is 
a doubly stretched function, one due to the continuously 
increasing relaxation time with time, and the other due 
to the spreading of the spectrum due to the heterogene- 
ity in domain sizes. The heterogeneity will cause the 
shape of the relaxation spectrum to change with temper- 
ature, as the apparent activation energy is proportional 
to the domain size. The shape will become flatter at 
lower temperatures. Viscoplastic data from the stress- 
strain curves of glassy polymers support the reduction 
of the power law exponent in the high frequency range 
as the temperature is lowered well below T^ [20]. 

8. Summary 

The Vogel equation can be derived from the tempera- 
ture dependence of the domain size for intermolecular 
cooperativity from the Adam-Gibbs formula in equi- 
librium. The measured entropy approaches zero at a 
temperature To that depends on the conformer size, 
though the molecules are still in the highly amorphous 
disordered state, i.e., high configurational entropy. In 
the cooperativity model, all conformers within a do- 
main are allowed to relax in unison. The internal degree 
of freedom within a domain is reduced to that of one 
conformer. As the measurable entropy is the configura- 
tional entropy per mole of conformers, it decreases 
more rapidly as the domain grows, than the conforma- 
tional entropy that describes disorder in polymer 
molecules. 



The heterogeneity of domain sizes are considered to 
be responsible for the distribution of relaxation times, as 
manifest in Kovacs' memory effect and Hodge's sub-T^g 
relaxation. A power law t~" is derived for the relaxation 
spectrum based on the free energy of domain size fluc- 
tuations. 

The van der Waals excess free volume is much greater 
than the Doolittle free volume. The latter is evaluated 
from the relaxation data, and its expansion coefficient. 
Off is not equal to A a , but consistently differs depending 
on the molecular size. From the size of the conformer, 
Tg can be estimated. For a single conformer to be able 
to relax free of the neighbors interference, it needs the 
free volume comparable to the occupied volume. It has 
been shown that the coefficient of thermal expansion for 
the free volume fraction is approximately equal to the 
ratio between the kT and A/x*, the latter being the en- 
ergy barrier for the bond rotation. The amount of free 
volume trapped at T^ is greater for the glass vitrified at 
a higher temperature, while its glassy modulus is lower. 
This is true for crosslinked polymers as well. 

A relationship has been derived relating the ongoing 
relaxation time to the aging time. The ratio of the char- 
acteristic relaxation time to the elapsed aging time is 
equal to the quantity kT/Afjiz, which is in the range of 
0.1 to 0.01 depending on the temperature and thermal 
history of the glassy state being aged. 
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